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Die folgenden Angaben sind den vom Anmelder eingereichten Unterlagen entnommen 

® Mehrdimensionale Monitore fiir statistische Prozesse 

© Es wird eine erweiterte Vorgehensweise mit partiellen 
kleinsten Fehlerquadraten (EPLS) fiir die Zustandsuber- 
wachung industrieller Prozesse beschrieben. Diese EPLS- 
Vorgehensweise stellt zwei statistische Uberwachungs- 
diagramme zur Verfugung, urn anomales Prozessverhal- 
ten festzustellen, und ebenfalls Beitragsdiagramme, urn 
dieses Verhalten zu diagnostizieren. Es wird eine theoreti- 
sche Analyse der EPLS-Uberwachungsdiagramme zur 
Verfiigung gestelrt, zusammen mit zwei Anwendungsstu- 
dien, urn zu zeigen, dass die EPLS-Vorgehensweise ent- 
weder empfindlicher ist oder eine einfachere Interpretati- 
on ermoglicht als herkommliche PLS. 
Verallgemeinerte Ergebnisse werden durch Konstruktion 
einervergroSerten Matrix folgender Form berechnet 

2»[YiX] # 

wobei X die Pradiktormatrix und Y die Responsmatrix ist 
und es wird eine Ergebnismatrix T n = T n - E n * konstruiert, 
bei welcher T n * und E n * folgende allgemeine Form auf- 
weisen: 

; t; = [Y;x]te:r.] l R„ 
| e; = (e.;fJ[bj;! 5 :3J 1 r ( , 

wobei die Spalten der Matrix T n * die verallgemeinerten 
t-Ergebnisse zur Verfugung stellen und die Spalten der 
Matrix E n * die verallgemeinerten restlichen Ergebnisse, 
wobei eine Einheitsmatrix MxM bezeichnet und 

PLS 

die PLS-Regressionsmatrix ist. 
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Beschieibung 

[0001] Die vorliegende Erfindung betriflt mehrdimensionale Monitore fur statistische Prozesse. Der BegrifF "Prozess" 
soli als allgemeiner BegrifF der Regeltheorie verstanden werdeo und allgemein gesteuerte bzw. geregelte Gerate, Anla- 
5 gen und derartige Systeme umfassen. 

1. EINFUHRUNG 
1.1 Hintergrund der Erfindung 

10 

[0002] Die Erf assung und Diagnose anomaler Situationen im Betrieb industrieller Prozesse stellt ein Problem mit be- 
trachtlicher Herausforderung dar, das sowohl im akademischen als auch industriellen Bereich grofie Beachtung findet. 
[Nimmo, 1995] erlauterte, dass allein die petrochemische Industrie in den USA bis zu 10 Milliarden Dollar pro Jahr ein- 
sparen konnte, wenn anomale Situationen erfasst, diagnostiziert und entsprechend behandelt werden konnten. Die Fol- 

15 gen der Tatsache, derartige Ereignisse nicht feststellen zu konnen, konnen von erhbhten Betriebskosten beim Ablauf ei- 
nes Prozesses bis zu Produktionsverlusten fuhren, infolge eines desastrosen Ausfalls der gesamten Anlage. 
[0003] Die Aufgabe der Erfassung und Diagnose industrieller Prozesse, kontinuierlich oder postenweise, ist schwierig. 
Dies liegt daran, dass bei industriellen Prozessen haufig eine groBe Anzahl an Prozessvariablen vorhanden sind, bei- 
spielsweise Temperaturen, Driicke, Russraten, Zusammensetzungen usw., die regelmafiig bis zu einigen tausend Malen 

20 pro Tag aufgezeichnet werden [Piovose, 1991] [Kosanovich, 1992]. Diese sehr grofie Datenmenge ist schwierig, einfach 
durch Beobachtung zu analysieren und zu interpretieren. Weiterhin ist es haufig der Fall, dass die Prozessvariablen hoch- 
korreliert sind [MacGregor, 1991] und daher die Anzahl an Freiheitsgraden innerhalb des Prozesses erheblich kleiner ist 
als die Anzahl beobachteter Prozessvariabler. Dies macht es selbst fur einen erfahrenen Bediener schwierig, die Wech- 
selwirkung von Ursache und Wirkung durch Beobachtung zu interpretieren. Allerdings sind in den aufgezeichneten Da- 

25 ten die Grundlagen zur Veranschaulichung des momentanen Zustands der Prozessoperation verborgen. Die Schwierig- 
keit besteht darin, diese Grundlagen aus den Daten abzuziehen. 

[0004] Um dieses Thema der Erfassung und Diagnose anzugehen, wurden Vorgehensweisen der mehrdimensionalen, 
statistischen Prozesssteuerung (MSPQ erfolgreich eingesetzt [Kresta, 1991], [MacGregor, 1995], [Kourti, 1995]. Die 
MSPC-Techniken sind darauf gerichtet, sukzessiv die Anzahl an Variablen zu verringern, die dazu erforderlich sind, si- 
30 gnifikante Variationen des Prozesses zu beschreiben. Die aufgezeichneten Daten werden daher in eine Gruppe von we- 
niger Variablen komprimiert, die daher besser handhabbar und interpretierbar sind. 

[0005] Eine derartige MSPC- Vorgehensweise sind die partiellen kleinsten Quadrate (PLS), die auf der Pionierarbeit 
von H. Wold Mitte der Sechziger beruhen [Gelade, 1988], Die ersten Veroffentlichungen von PLS wurden 1966 vorge- 
stellt [Wold, 1966a; 1966b]. Das PLS-Verfahren identifiziert eine parametrische Regressionsmatrix auf der Grundlage 

35 von Pradiktor- und Respons-Matrizen, die aus Referenzdaten des Prozesses aufgebaut werden. Die Pradiktormatrix be- 
steht aus den Signalen der manipulierten und gemessenen Storungs- oder Ursachen variablen des Prozesses (Pradiktorva- 
riable), wogegen die Responsmatrix aus dem Gesteuerten oder Ergebnisvariablen des Prozesses (Responsvariablen) be- 
steht. DerPLS-Algorithmus zerlegtdie Pradiktormatrizen und Responsmatrizen in Komponentenmatrizen des Ranges 1. 
Jede Komponentenmatrix besteht aus einem Vektorprodukt, bei welchem ein Vektor die Variation (Ergebnisvektor) und 

40 der andere den Beitrag (Belastungsvektor) des Ergebnisvektors entweder zur Pradiktormatrix oder zur Responsmatrix 
beschreibt. Die Zeriegung ist eine Iteration, bei welcher ein Paar von Komponentenmatrizen (eine fur den Pradiktor und 
eine fur die Respons) in jedem Iterationsschritt berechnet wird. Die Regressionsmatrix wird bei jedem Iterationsschritt 
als Ergebnis dieser Zeriegung aktualisiert. Die Datenreduktion wird dadurch erzielt, dass die \fcriation der Pradiktor- und 
Responsvariablen bis zur kleinsten Anzahl an Ergebnisvektoren herunterkomprirniert wird, die wirksam das Prozess ver- 

45 halten beschreiben konnen. Die Auswahl der Anzahl an Komponentenmatrizen, die beibehalten werden mussen, stellt ei- 
nen Kompromiss zwischen dem Maximieren der Variation, die in den Pradiktor- und Responsmatrizen erlautert wird, 
und dem Minimieren der Anzahl an Komponentenmatrizen dar. Kreuzuberprufung [Wold, 1978] wird am haufigsten 
dazu verwendet, die Anzahl zuriickzuhaltender Komponentenmatrizen zu definieren, vgl. beispielsweise [MacGregor, 
1991, 1995], [Morud, 1996]. 

50 [0006] [MacGregor, 1995] und [Wise, 1996] gaben an, dass die PLS-Zerlegung der Pradiktormatrix fur die Zustands- 
uberwachung kontinuierlicher, industrieller Prozesse verwendet werden kann. Weiterhin betonen sie, dass diese Zerie- 
gung ahnlich einer Hauptkomponentenanalyse (PCA) der Pradiktormatrix ist. Die PLS-Zerlegung der Pradiktormatrix 
gestattet die Berechnung von zwei Statistiken. Die erste Statistik (die T-Quadratstatistik) beschreibt die Variation der 
Pradiktormatrix, die signifikant zur Vorhersage der Responsvariablen ist. Im Gegensatz hierzu entspricht die zweite Sta- 

55 tistik (die Q-Statistik) der Variation bei der Pradiktormatrix, die insignifikant zur \fcrhersage der Responsvariablen ist. 
Beide Statistiken konnen in statistischen Uberwachungsdiagrammen in Abhangigkeit von der Zeit aufgetragen werden. 
Diese Vorgehensweise wird nachstehend als "Vorgehensweise I" bezeichnet 

[0007] Eine weitere Vorgehensweise zur Nutzung von PLS als ein Zustandsuberwachungswerkzeug wird beispiels- 
weise in [MacGregor, 1991] und [Kresta, 1991] diskutiert. Bei dieser \forgehensweise werden verschiedene statistische 
60 Plots dazu verwendet, ein anomales Prozessverhalten festzustellen und zu diagnostizieren. Diese Plots sind: 

- x-y-Plots des quadratierten \forhersagefehlers der Responsvariablen in Abhangigkeit von den Ergebniswerten je- 
des Ergebnisvektors, der die Pradiktorvariablen reprasentiert (Uberwachungsdiagramme), 

- Rots jeder Kombination der beiden Ergebnisvektoren, welche die Pradiktorvariablen reprasentieren (Streuungs- 
65 Plots), und 

- Plots des quadrierten Vorhersagefehlers, entweder der Pradiktor- oder Responsvariablen, in Abhangigkeit von 
der Zeit (SPE-Diagramme). 
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[0008] Diese Vorgehensweise wind in der nachstehenden Beschreibung als "Vorgehensweise 11" bezeichnet. 

1.2 Zusammenfassung derErfindung 

[0009] In der folgenden Beschreibung wird eine Ausweitung des Standard-PLS-Algorithmus, nachstehend als "erwei- 5 
terter FLS" oder "EPLS" bezeichnet, fur kontinuierliche Prozesse eriautert. Diese Erweiterung fuhrt zu der Bestimmung 
zweier neuer PLS-Ergebnisse auf der Grundlage der Ergebnisvektoren der Pradiktormatrix. Die neuen PLS-Ergebnis- 
vektoren werden als verallgemeinerte Ergebnisvektoren bezeichneL Der erste verallgemeinerte Ergebnisvektor be- 
schreibt eine signifikante Variation des Prozesses einschlieBlich der Pradiktor- und Responsvariablen. Der zweite verall- 
gemeinerte Ergebnisvektor reprasentiert den Vorhersagefehler des PLS-Modells und Residuen der Pradiktormatrix. Die 10 
EPLS- Vorgehensweise fuhrt zu Uberwachungsdiagrammen fiir T-Quadrat und Q, die ahnlich jenen sind, die bei PCA er- 
halten werden, wenn sowohl Pradiktor- als auch Responsvariablen mittels PCA analysiert werden. Dies unterscheidet 
sich von der ublichen PIS- Vorgehensweise, welche nur die Pradiktorvariablen analysiert und daher keinen Aufschluss 
uber das Verhalten der Responsevariablen ermoglicht, es sei denn, in dem Prozess ware eine Ruckkopplung vorhanden. 
[0010] Der Vorteil der EPLS-Cberwachungsdiagramme besteht daher darin, dass sie eine Variation der Pradiktor- und 15 
Responsvariablen zusammen mit ihren Residuen reprasenueren. Dies verbessert die Oberwachungsdiagramme der \br- 
gehensweise I, welche nur die Variation und Residuen der Pradiktorvariablen beschreiben. Im Gegensatz zur \forgehens- 
weise II, stellt EPLS die Fahigkeit zur Verfugung, den Prozess auf nur zwei Diagrammen zu uberwachen, im Gegensatz 
zu bisher, als die Anzahl an Diagrammen von der Anzahl an Komponentenmatrizen abhing. 

[0011] GemaB einer Zielrichtung der vorliegenden Erfindung umfasst ein Verfahren fur die Konstruktion und das Kon- 20 
figurieren eines mehrdimensionalen Monitors fur statistische Prozesse durch eine Vorgehensweise mit partiellen klein- 
sten Quadraten die Konstruktion der Pradiktor- und Responsmatrizen des Prozesses aus Referenzdaten, wobei die Pra- 
diktormatrix aus Signalen der manipulierten und gemessenen Storungs- oder Ursachenvariablen des Prozesses (Pradik- 
torvariablen) besteht, und die Responsmatrix aus den gesteuerten oder Auswirkungsvariablen des Prozesses (Respons- 
variablen) besteht, Zerlegen der Pradiktor- und Responsmatrizen in Komponentenmatrizen des Rangs 1, wobei jede der 25 
Komponentenmatrizen aus einem Vektorprodukt besteht, bei welchem ein Vektor (der Ergebnisvektor) die Variation be- 
schreibt, und der andere (der Belastungsvektor) den Beitrag des Ergebnisvektors zu der Pradiktor- oder Responsmatrix 
beschreibt, und die Zerlegung durch die Erzeugung einer parametrischen Regressionsmatrix durchgefuhrt wird, auf der 
Grundlage von Iterationen der Zerlegung der Pradiktor- und Responsmatrizen, und ist gekennzeichnet durch die Erzeu- 
gung eines ersten verallgemeinerten Ergebnisvektors, der jede signifikante Variation des Prozesses beschreibt, ein- 30 
schlieBlich Variationen der Pradiktor- und Responsvariablen, sowie eines zweiten verallgemeinerten Ergebnisvektors, 
welcher den Vorhersagefehler des Modells der partiellen kleinsten Quadrate und Residuen der Pradiktormatrix reprasen- 
tiert. 

[0012] Vorzugsweise werden die verallgemeinerten Ergebnisse durch Konstruktion einer vergrbBerten Matrix berech- 
net, die hiermit mit Z bezeichnet ist und folgende Form aufweist: 35 

z=[y:x], 

wobei X die Pradiktormatrix und Y die Responsmatrix ist, und durch Konstruktion einer Ergebnismatrix T„ = T n *-En*, 
wobei T n * und E n * folgende allgemeine Form aufweist: 40 

t; = [y;x][b&:3)X 

wobei die Spalten der Matrix T n * die verallgemeinerten t-Ergebnisse angeben, und die Spalten der Matrix En* die ver- 
allgemeinerten, restlichen Ergebnisse, wobei 3 eine Einheitsmatrix M x M bezeichnet, und B}^ die PLS-Regressions- 
matrix ist 

[0013] GemaB einer zweiten Zielrichtung der Erfindung stellen wir einen mehrdimensionalen Monitor fiir statistische 50 
Prozesse zur Verfugung, der gemaB der erste Zielrichtung der Erfindung konstruiert und konfiguriert ist, und der so aus- 
gebildet ist, dass er anomales Prozessverhalten durch Analyse der Residuen der Responsvariablen identifiziert. 
[0014] GemaB einer dritten Zielrichtung der Erfindung stellen wir ein Verfahren zur Oberwachung eines Prozesses zur 
Verfugung, welches das Konfigurieren eines mehrdimensionalen Monitors fur statistische Prozesse mit dem Verfahren 
gemaB der ersten Zielrichtung der Erfindung umfasst, und zum Identifizieren eines anomalen Prozessverhaltens, zumin- 55 
dest teilweise durch Analyse der Residuen der Responsvariablen. 

[0015] Die Erfindung wird nachstehend, nur beispielhaft, unter Bezugnahme auf die beigefugten Figuren beschrieben. 
Es zeigt: 

[0016] Fig, 1 eine schematische Darstellung der katalytischen Huidkrackeinheit, 

[0017] Fig. 2 eine schematische Darstellung eines FlieBbettreaktors und seiner benachbarten EinheiL 60 
[0018] Figuren fur die katalytische Huidkrackeinheit ohne Regelungsriickkopplung in der Pradiktormatrix: 
[0019] Fig, 3 Statistikuberwachungsdiagramme fur normale Betriebsdaten (die oberen Diagramme geben die FLS- 
Uberwachungsdiagramme an, PLS-T - und Q-Statistiken, und die unteren Diagramme zeigen die EPLS-Uberwachungs- 
diagramme, EPLS-T 2 - und Q-Statistiken), 

[0020] Fig. 4 statistische Oberwachungsdiagramme fur die ungemessene Stoning (die oberen Diagramme geben die 65 
PLS-t)berwachungsdiagramme an, PLS-T 2 - und Q-Statistiken, und die unteren Diagramme zeigen die EFLS-Oberwa- 
chungsdiagramme, EPLS-T 2 - und Q-Statistiken, 

[0021] Fig. 5 Fehlerbeitragsdiagramm fur den Zeitpunkt 11460 Minuten. Die Or und die COKonzentration in dem 
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Dichtgasfluss weisen den groBten Vorhersagefehler auf, 

[0022] Fig. 6 statistische Uberwachungsdiagramme fur die Anderung des Flusses des regenerierten Katalysators in den 
Reaktor (die oberen Diagramme zeigen die PLS-t)berwachungsdiagramme, PLS-? 2 - und Q-Statistiken, und die unteren 
Diagramme zeigen die EPLS-t)berwachungsdia§ramme, EPLS-T 2 - und Q-Statistiken), 
5 [0023] Fig. 7 Fehlerbeitragsdiagramm fur die Anderung des Flusses des regenerierten Katalysators zum Reaktor zum 
Zeitpunkt 19357 Minuten. Der Steigrohr-Katalysator-Pegel und die (VKonzentration in dem Gichtgas werden am stark- 
sten beeinflusst 

[0024] Figuren fur die katalytische Fluidkrackeinheit mit Regelungsriickkopplung in der Pradiktormatrix: 
[0025] Fig. 8 statistische t)berwachungsdiagramme fur die ungemessene Stoning (Verkokungsfaktor) ; Pradiktorvaria- 
10 blen umfassen das Feuchtgaskompressor-Saugventil (obere Diagramme sind die PIi>-t)berwachimgsaiagramme, PLS- 
T 2 - und Q-Statistiken, und die unteren Diagramme sind die EPIi^Uberwachungsdiagrarrinie, EPLS-T 2 - und Q-Statisti- 
ken), zum Zeitpunkt 1 9357 Minuten. 

[0026] Fig. 9 statistische Uberwachungsdiagramme fur die Anderung des Flusses des regenerierten Katalysators zum 
Reaktor; Pradiktorvariablen umfassen das Feuchtgaskompressor-Saugventil (obere Diagramme sind die PLS-Uberwa- 
15 chungsdiagramme, PLS-T 2 - und Q-Statistiken, und die unteren Diagramme sind die EPI^-t)berwachungsdiagrarnme, 
EPLS-T 2 - und Q-Statistiken). 
[0027] Figuren des FlieBbettreaktionsprozesses: 

[0028] Fig. 10 Statistikuberwachungsdiagramme fur Normalbetriebsdaten (oberes Diagramm ist das PLS-Uberwa- 
chungsdiagramm, PLS-T 2 - und Q-Statistiken, und die unteren Diagramme sind die EPLS-Uberwachungsdiagramme, 
20 EPLS-T 2 - und Q-Statistiken), 

[0029] Fig. 11 statistische Oberwachungsdiagrarrirne fur die ungemessene Stoning (oberes Diagramm ist das PLS- 
Ubemachungsdiagramm, PLS-T 2 - und Q-Statistiken, und die unteren Diagramme sind die EPLS-Uberwachungsdia- 
gramme, EPLS-T 2 - und Q-Statistik), 

[0030] Fig. 12 EC-Diagramme fur den Dampfdmckumschlag zu den Zeitpunkten 1500 Minuten (oberer linker Plot), 
25 1501 Minuten (unterer linker Plot) und 1502 Minuten (oberer reenter Plot), 

[0031] Fig. 1 3 statistische Dberwachungsdiagranime fur ein anomales Verhalten eines der Rohre (oberes Diagramm ist 
das PLS-t)berwachungsdiagramm, PLS-T 2 - und Q-Statistik, untere Diagramme sind die EPLS-Oberwachungsdia- 
gramme, EPLS-T 2 - und Q-Statistiken), 

[0032] Fig. 14 EC-Diagramme fur das Fluidisierungsproblem in einem der Rohre zu den Zeitpunkten 436 Minuten 
30 (oberes linkes Diagramm), 888 Minuten (oberes rechtes Diagramm), 905 Minuten (unten links) und 910 Minuten (unte- 
res rechtes Diagramm). 

[0033] Um die Niitzlichkeit der EPLS-t)berwachungsdiagramme zu demonstrieren und einen Vergleich mit der Vorge- 
hensweise I durchzufuhren, werden nunmehr zwei Fallbeispiele betrachtet. Fiir jedes werden zwei typische Bedingungen 
fur anomales Verhalten erzeugt, welche die Auswirkung einer ungemessenen Stoning und ebenso eine "interne" Ande- 

35 rung des Prozessverhaltens beschreibt. Die Prozessbeispiele sind die Simulation einer katalytischen Fluidkrackeinheit 
(FCCU), eingefuhrt von [McFarlane, 1993], und ein realer industriellerProzess, der zwei unterschiedliche Losungsmittel 
als Ergebnis einer komplizierten chemischen Reaktion erzeugt, die in einem FlieBbettreaktor durchgefuhrt wind. 
[0034] Die folgende Beschreibung ist folgendermaBen organisiert. Im nachstehenden Abschnitt 2 werden sowohl der 
Standard-PLS-Algorithmus als auch der neue EPLS-Algorithmus beschrieben und verglichen. Im Abschnitt 3 werden 

40 die Zustandsiiberwachungsstatistiken eingefuhrt, welche diesen beiden \fcrgehensweisen zugeordnet sind. Abschnitt 4 
gibt die Anwendungsfalluntersuchungen an, um Beispiele fur die \forteile von EPLS anzugeben. 

2. ALGORITHMEN PARTELLER KLEINSTER QUADRATE 

45 2. 1 Der Standard-PLS-Algorithmus 

[0035] Die S tandard-PLS-Identifizierungstechnik beruht auf der Zerlegung der Pradiktormatrix Xb e R 1 ^ und der Re- 
sponsmatrix YqER 1004 in eine Summe von KomponentenmaUizen des Rangs 1 [Gelade, 1986]. Beide Matrizen enthal- 
ten K Datenpunkte, die Pradiktormatrix besteht aus M Variabien, und die Responsmatrix aus N Variablen. Beide Matri- 
50 zen werden gewohnlich vor der Identi fizierungsprozedur mittenzentriert und geeignet skaliert. Die Zerlegung beider Ma- 
trizen erfolgt folgendermaBen: 

55 . (1) 

60 wobei X und Y die Komponentenmatrix der Pradiktor- bzw. Responsmatrix ist GemaB Gleichung (1) kbnnen die Matri- 
zen des Rangs Eins als Vektorprodukt zwischen ti und % definiert als Ergebnisvektoren oder latente Variablen (LVs), und 
von pi und oj, definiert als Belastungsvektoren, berechnet werden. M ist gleich der Anzahl an Pradiktorvariablen, und ^ 
reprasentiert den Vorhersagefehler des Prozessmodells. Es wird darauf hingewiesen, dass dann, wenn alle Komponen- 
tenmatrizen eingeschlossen sind, die Pradiktonnatrix gleich den Matrixzerlegungen ist. V\fenn nur n Komponentenmatri- 

65 zen eingeschlossen sind, dann wird Gleichung (1) zu: 
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(2J 

Y 0 = Z Y, + E„ = £ u.qr + E„ = U. • Ql + E„ 

#»l J-l 

wobei F n die Residuen der Pradiktormatrix bezeichnen. Die vorhergesagten u-Ergebnisse U| konnen durch die folgende 
Multiplikation bestimmt werden: 10 

U n = [t l b l :t D b n ]=T n diag{b n }, (3) 

wobei diag{b n } eine Diagonalmatrix ist, welche die Regressionskoeffizienten b\ des Ergebnismodells aufeinanderfol- 
gend enthalL Eine theoretische Untersuchung des PLS-Algorithmus findet sich im Anhang 1. Unterschiedliche \forge- 15 
hensweisen wurden eingefuhrt, um die Ergebnis- und Belastungsvektoren zu bestimmen, und dies sind der LSQR-Algo- 
rithmus [Manne, 1987], der NIPALS-Algorithmus [Geladi, 1986], der SIMPLS-Algorithmus [de Jong, 1993] und andere. 
[0036] Infolge der Tatsache, dass industrielle Prozesse haufig stark korrelierte Prozessvariablen aufweisen, konnen nur 
einige wenige LVs dazu erforderlich sein, den Hauptanteil der Prozess variation zu beschreiben. Im Gegensatz hierzu be- 
treffen die verbleibenden Paare der LVs im Wesentlichen Rauschen und nicht-signifikante Variationen bei Xo und Y 0 20 
[Gelade, 1986], [MacGregor, 1991] und [Wise, 1996]. Um die Anzahl zuriickzuhaltender LVs zu bestimmen, wurden 
Kreuzuberpriifung [Wold, 1978] und die Analyse der Varianz (ANOVA) [Jackson, 1991] diskutiert. 



2.2 Ableitung des EPLS-Algorithmus und der verallgemeinerten Ergebnisse 

25 

[0037] Der EPLS-Algorithmus erzeugt Ergebnisse, welche die Variation der Pradiktor- und Responsvariablen repra- 
sentieren, sowie deren Residuen, und diese werden als verallgemeinerte Ergebnisse bezeichnet. Diese Ergebnisse bilden 
die Grundlage fur eine wirksamere Prozesszustandsuberwachung als bei den vorhandenen \forgehensweisen, die haupt- 
sachlich auf Ergebnissen beruhen, welche nur die Variation der Pradiktorvariablen beschreiben. Die verallgemeinerten 
Ergebnisse werden berechnet, nachdem die Gewichtungs- und Belastungsmatrizen bestimmt wurden (siehe Anhang 1), 30 
und beruhen auf der VergroBerung der Responsmatrix zu der Pradiktormatrix. Die vergroBerte Matrix wird mit Z be- 
zeichnet und ist folgendermaBen definiert: 

z=[y;x] (4). 

35 

[0038] Es wird darauf hingewiesen, dass der Index 0 bei beiden Matrizen weggelassen isL Dies liegt daran, dass die 
Ableitung der verallgemeinerten Ergebnisse auf dem Standard-PLS-Algorithmus beruht und der Druckablassvorgang 
nicht erneut durchgefuhrt werden muss. 

[0039] Aus der Gleichung (4) erhalt man durch Subtraktion der vorhergesagten Responsmatrix Y n und der rekonstru- 
ierten Pradiktormatrix X n , wobei eine Anzahl von n an LVs beibehalten wird, den folgenden Ausdruck: 40 

[Y:X]-|_Y n iX n J=CE a :F 0 ], (5) 

wobei En der Vorhersagefehler der Responsmatrix ist, und F n die Residuen der Pradiktormatrix bezeichnen. Durch Ein- 
setzen der Gleichungen (2) und (3) kann Gleichung (5) folgendermaBen umgeschrieben werden: 45 

[Y:X]-T J ,l^{b„}Ql:P;J S! [EjF.]. (6) 



[0040] Wie im Anhang 2 gezeigt, ist das Matrixprodukt P T n mit der PLS-Regressionsmatrix B^, wobei n LVs zuriick- 50 
gehalten werden, gleich dem Matrixprodukt der Diagonalmatrix diag{bn} und Qq T . Das Integrieren dieses Ergebnisses in 
Gleichung (6) ergibt: 



[YiX]-T„P;lB« P I»;3j=[EjFj. (7) 

wobei 3 eine Einheitsmatrix M x M bezeichnet. Die Ausfuhrung einer Multiplikation von Gleichung (7) mit dem verall- 
gemeinerten Inversehen von [B PL s : 3] ergibt: 

[Y:X][B{& :3f - T,P; = [E„ iF.fe itf , {g) 

wobei J- das verallgemeinerte Inverse bezeichnet Wie im Anhang 3 gezeigt, fuhrt die Multiplikation von Gleichung (8) 
mit R n (siehe Anhang 1) zu einer Formel zur Berechnung der Ergebnisse der Pradiktormatrix T n : 

T. = [YJxiB&iofc - [E.:F.Jte:3] l R < , . (9) 

[0041] In Gleichung (9) ist die Ergebnismatrix T n gleich der Differenz von zwei Matrizen. Die erste Matrix bezieht 
sich auf die Pradiktor- und die Responsmatrix, und die zweite Matrix hangt von dem Vbrhersagefehler der Responsma- 
trix und den Residuen der Pradiktormatrix ab. Die Matrix [En i FJ wird als die vergroBere Residuenmatrix von F n be- 
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zeichnet. 

[0042] Definiert man die Matrix 

als CS*s, die sich ergebende Matrix [Y ! X]Cj^ als T n *, und das Matrixprodukt [En i FJC^ als En*, so vereinfacht sich 
Gleichung (9) folgendermaBen: 

T n E T n *-En*. (10) 

10 

[0043] Die Spalten der Matrix T n * werden als die verallgemeinerten t-Ergebnisse bezeichnet, wahrend die Spalten der 
Matrix En* als verallgemeinerte Residuenergebnisse bezeichnet werden. Fur die Prozesszustandsuberwachung stellt 
Gleichung (10) Ergebnisse zur Verfiigung, welche die Prozessvariation beschreiben, die innerhalb der Absolutwerte der 
Pradiktor- und der Responsmatrix enthalten sind, sowie die Vorhersagefehlermatrix E„, als auch die Residuen der Pra- 
15 diktormatrix F n . 

[0044] Der nachste Abschnitt beschreibt die Ableitung der Statistiken fur T„* und E n *, die in Abhangigkeit von der 
Zeit bei eindimensionalen Uberwachungsdiagrammen aufgetragen werden konnen. 

2.3 Vergleich mit den Vorgehensweisen I und II fur die Prozesszustandsuberwachung 

20 

[0045] Wie voranstehend erwahnt, beruhen die vorhandenen Vorgehensweisen zur Prozesszustandsuberwachung 
hauptsachlich auf den t-Ergebnissen. [MacGregor, 1995] und [Wise, 1996) fuhrten aus, dass bei der Vorgehensweise I die 
weggefallenen und die beibehaltenen t-Ergebnisse die Grundlage fur zwei Uberwachungsdiagramme bilden, die in dem 
nachsten Unterabschnitt diskutiert werden. Obwohl erfolgreiche Anwendungen der \fcrgehensweise I diskutiert wurden, 
25 beispielsweise [Kourti, 1995], [Wise, 1996] und [Morud, 1996], stellen sie nicht notwendigerweise jede Art eines an- 
omalen Prozessverhaltens fest. Dies gilt insbesondere, wenn: 

1. anomales Prozessverhalten hauptsachlich die Responsvariablen beeinflusst, die nicht durch eine geschlossene 
Regelschleife geregeit werden. In diesem Fall breitet sich das anomale Verhalten nicht bis zu den Pradiktorvariablen 

30 mittels Regelungsriickkopplung aus, und bleibt daher unentdeckt. Mit EPLS wird die Variation der Responsvaria- 

blen deutlich. 

2. Die Responsvariablen sind stark korreliert, jedoch die Pradiktorvariablen nicht. In diesem Fall kann nur ein sta- 
tistisches Diagramm fiir die Standard-PLs- Vorgehensweise erhalten werden, das Diagramm fur quadriertes T. Bei 
EPLS bleiben beide Diagramme relevant, unabhangig von der Anzahl beibehaltener LVs. 

35 

[0046] Die Vorgehensweise II verlasst sich auf Streuungsplots und x-y-Diagramme von SPE in Abhangigkeit von ver- 
schiedenen t-Ergebnissen und SPE-Diagrammen, beispielsweise [MaxGregor, 1991] und [Kresta, 1991]. Wenn derPro- 
zess aus einer groBen Anzahl stark korrelierter Prozessvariablen besteht, beispielsweise einhundert oder mehr, kann die 
Anzahl erforderlicher Ergebnisse allerdings immer noch groB sein, damit signifikante Prozessvariationen erf asst werden. 
40 Daher kann bei diesen anderen Vorgehensweisen eine groBe Anzahl an Diagrammen erforderlich sein, und wird die Si- 
tuation mtihsam zu analysieren. Im Gegensatz hierzu sind bei EPLS nur zwei Diagramme erforderlich, unabhangig von 
den Dirnensionen des Problems und von der Anzahl ausgewahlter LVs. 

[0047] Man konnte auch das SPE-Diagramm von Responsvariablen zusatzlich zur Vorgehensweise I dazu verwenden, 
die voranstehenden Nachteile zu iiberwinden. Hierbei wurden jedoch zumindestdrei Ubemachungsdiagramme benotigt, 
45 und wird die Variation der Responsvariablen nicht in einem dieser Diagramme akkumuliert Im Gegensatz hierzu erfor- 
dern die verallgemeinerten Ergebnisse nur zwei Uberwachungsdiagramme, und erfasst eines der generalisierten Ergeb- 
nisse die Variation der Responsvariablen. SchlieBlich wird darauf hingewiesen, dass die Verwendung der generalisierten 
Ergebnisse zur Prozesszustandsuberwachung ahnlich der Art und Weise ist, in welcher PCA zur Uberwachung indu- 
strieller Prozesse verwendet wird [Jackson, 1991]. 

50 

3. STATISTIK DER PLS- UND EPLS- VORGEHENSWEISE 
3.1 Statistik fur PLS 

55 [0048] Bei der Vorgehensweise I konnen zwei statistische Uberwachungsdiagramme erhalten werden, auf der Grund- 
lage der Zerlegung der Pradiktormatrix. Das erste Uberwachungsdiagramm steht in Beziehung zu den beibehaltenen t- 
Ergebnissen und beschreibt einen signifikanten Beitrag zur \forhersage der Responsmatrix. Das zweite Diagramm ist der 
Variation der Pradiktormatrix zugeordnet, die von den weggefallenen t-Ergebnissen erfasst wird. Die weggefallenen t- 
Ergebnisse beschreiben nicht-signifikante und unkorrelierte Beitrage zur \forhersage der Responsmatrix. Allerdings 

60 muss in einem Fall, in welchem samtliche Pradiktorvariablen significant zur Variation der Responsvariablen beitragen, 
jedes t-Ergebnis beibehalten werden. Daher bleiben keine t-Ergebnisse zur Berechnung des zweiten Oberwachungsdia- 
gramms ubrig. 

[0049] Das erste Uberwachungsdiagramm beruht auf einer Statistik, die als PLS-T^Statistik bezeichnet wird, und das 
zweite Dberwachungsdiagramm betrifft eine Statistik, die als PLS-Q-Statistik bezeichnet wird. Diese beiden Statistiken 
65 sind folgendermaBen definiert 
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(11) 



wobei plsT\ und PLsQk die PLS-T 2 - bzw. die Q-Statistik bezeichnen. Weiterhin bezeichnet t^ den Wert des i-ten t-Er- 10 
gebnisses zum Zeitpunkt k, und bezeichnet die Standardabweichung des i-ten Ergebnisvektors der Pradiktorvariablen 
der Bezugsdaten. f^i reprasentiert die Residuen der i-ten Pradiktorvariablen zum Zeitpunkt k, und x Ci ist die Standardab- 
weichung der i-ten Residuenvariablen der Bezugsdaten. Die Bezeichnungen T 2 und Q wurden entsprechend der T 2 - und 
Q-Statistik von Hotelling ausgewahlt, die bei PCA verwendet werden. Jede Statistik kann in einem Oberwachungsdia- 
gramm in Abhangigkeit von der Zeit aufgetragen werden. Es wird darauf hingewiesen, dass die Normierung von ^ und 15 
fki unbedingt erforderlich ist, um eine empfindliche Statistik zur Verfugung zu stellen. Falls dies nicht erfolgt, dominieren 
die t-Ergebnisse mit einer groBen Varianz, normalerweise die wenigen ersten, den Ergebniswert von PLS-T 2 , und die Re- 
siduen der Pradiktorvariablen, die eine groBe Varianz aufweisen, verdecken Residuen mit relativ kleiner Varianz. Ein 
Fehlerzustand, der hauptsachlich die t-Ergebnisse oder Residuen beeinflusst, welche kleine Anderungen aufweisen, kann 
in diesem Fall unentdeckt bleiben. Daruber hinaus kann in Bezug auf die Summe stochastischer \fcriablen mit Mittelwert 20 
Null und Einheitsvarianz (Chi-Quadrat-Verteilung) die statistische Bestimmung von Schwellenwerten verwendet wer- 
den. 

[0050] Wenn auBergewohnlich groBe PLS-T 2 - Werte auftreten, dann ist die gesamte Prozessvariation ungewohnlich 
groB, verglichen mit den Bezugsdaten des Prozesses. Dies bedeutet, dass das allgemeine Prozessverhalten sich betracht- 
tich geandert hat, oder dass sich der Prozess zu einem neuen Betriebsbereich bewegt hat. Im Gegensatz hierzu zeigen un- 25 
gewohnlich groBe PLS-Q-Werte an, dass sich die Beziehungen zwischen den Pradiktorvariablen in Bezug auf die Bezie- 
hung geandert haben, die bei den Bezugsdaten vorherrscht. 



3.2 Statistik fur EPLS 

[0051] Im Vergleich zu den Standard- t-Ergebnissen sind die statistischen Eigenschaften der verallgemeinerten Ergeb- 
nisse nachstehend zusammengefasst. 



[0052] Um die verallgemeinerten Ergebnisse fur die Prozesszustandsuberwachung zu analysieren, ist es wiinschens- 
wert, statistisch unabhangige Ergebnisse zu haben, was Orthogonalitat erfordert Um orthogonale Ergebnisse zu erzielen, 
kann eine Singularwertzerlegung (SVD) [Golub, 1996] der verallgemeinerten Ergebnisse eingesetzt werden, woraus sich 
ergibt: 

X* - v lT " ) A |T " ) \V tT ' )7 

" " " . (12) 



30 



1. Die verallgemeinerten t-Ergebnisse und ebenso die verallgemeinerten restlichen Ergebnisse sind um den Mittel- 
wert zentriert, wenn die Spalten in der Pradiktor- und der Responsmatrix in Bezug auf den Mittelwert vor der PLS- 35 
Identifizierung zentriert wurden. 

2. Die t-Ergebnisse des Standard-PLS-Algorithmus sind orthogonal [Hoskuldsson, 1988]. Im Gegensatz zu den 
Standard-t-Ergebnissen sind bei EPLS beide Ergebnistypen nicht orthogonal, unabhangig von der Anzahl zuriick- 
gehaltener LVs. Ein Beweis findet sich im Anhang 4. 



40 



45 



50 



wobei Vn^, A n ^ und W n ^ die SVD der verallgemeinerten t-Ergebnisse beschreiben und V.^, und W B ^ 
die SVD der verallgemeinerten, restlichen Ergebnisse reprasentierten. Die Dimensionen dieser Matrizen sind wie folgt: 
und V n (E *> sind Kxn-Matrizen, und A n <r*>, A n ^\ W/ 1 ^ und sind nxn-Matrizen. Die Spalten der Matrizen 
Vn 0 '^, V n (E *>, und W B ^ sind orthogonal, und Af*** und sind vom Diagonaltyp. Auf der Grundlage der 

SVD der verallgemeinerten Ergebnismatrizen der Bezugsdaten stellt die folgende Beziehung orthogonale Ergebnisse zur 55 
Verfugung: 

vi r) =T;wf'Ai r '"ViM 

V^n =E:w ,E-, A ^ ) -'^- T • (13) 

wobei Vn* 7 ** und V^** orthogonale T n *- und E^-Ergebnisse mit Einheitsvarianz reprasentieren. EinschlieBlich Glei- 
chung (13) konnen die orthogonalen Ergebnisse V^*' und Vn^ direkt aus den vergroBerten Daten und der Fehlerma- 
trix, [Y i X] und [En : FJ, folgendermaBen berechnet werden: 

65 
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[0053] Fiir die verallgemeinerten Ergebnisvektoren des i-ten Datenpunktes, i v (T * ) und iv,,^ kann die Summe der qua- 
drierten Elemente dazu verwendet werden, eine eindimensionale Statistik fiir jeden Vektor zu definieren. Diese Statisti- 
ken werden als t* • T 2 und E « • T 2 bezeichnet. t» • T 2 und E * • T 2 reprasentieren die EPLS-T 2 - und EPLS-Q-Statistik und 
10 sind folgendermaBen definiert 



T 2 


= V v« r > 











(15) 



[0054] Unter der Annahme, dass yV* 1 "** und y v (E * J stochastische Variablen sind, weisen beide Statistiken eine Chi-Qua- 
20 drat-Verteilung mit n Freiheitsgraden auf, welche die Vertrauensgrenzen zur Untersuchung dazu zur Verfugung stellt, ob 
der Prozess sich normal oder anomal verhalt. Die Vertrauensgrenzen sind normalerweise so ausgewahlt, dass sie 95% 
und 99% der Population enthalten (EPLS-T 2 - oder Q-Werte). Wenn ein neuer EPLS-T 2 - oder Q-Wert unterhalb der 
Grenze liegt, so wird die Hypothese akzeptiert, dass sich der Prozess normal verhalt, anderenfalls wind sie zuriickgewie- 
sen, und besteht die akzeptierte Hypothese darin, dass sich der Prozess anomal verhalt. 
25 [0055] Anomal groBe EPLS-T* 2*- und/oder Q-Werte konnen auftreten, wenn die Beziehung zwischen den Pradiktor- 
und den Responsvariablen, reprasentiert durch die parametrische Regressionsmatrix, sich geandert hat (beispielsweise 
ein zeitabhangiger Prozess), oder sich die Storungsstatistiken geandert haben. Andere Griinde konnen darin bestehen, 
dass der Prozess an einem unterschiedlichen Arbeitspunkt arbeitet, eine exzessive \fcriation des Prozesses aufgetreten ist, 
die nicht in den Bezugsdaten vorhanden war, oder dass ein anomales Prozessverhalten aufgetreten ist. Der Hypothesen- 
30 test ist daher ein Vergleich zwischen dem momentanen Prozessbetrieb und dem Prozessbetrieb, der durch die Bezugsda- 
ten wiedergegeben wird. Es wird darauf hingewiesen, dass die Bezugsdaten den Prozess im Normalbetrieb beschreiben, 
und jede Variation wiedergeben miissen, die im Normalbetrieb auftreten kann, da anderenfalls der statistische Hypothe- 
sentest zu empfindlich wird. 

[0056] [Wise, 1996] betonte, dass insbesondere die T 2 -Statistik nicht normal verteilt sein kann. [Dunia, 1996] unter- 

35 suchte den Einfluss eines exponentiell gewichteten, beweglichen Mittelwerts (EWMA) auf die Q-Statistik, welche PCA 
verwendet. Es stellte sich heraus, dass die mittlere Laufzeit (ARL), namlich die mittiere Zeit, die verstreicht, bis ein an- 
omales Prozessverhalten festgestellt wird, zur Feststellung von Fehlerzustanden bei Sensoren dadurch verringert werden 
kann, dass die EWMA-Q-Statistik verwendet wurde. In der vorliegenden Beschreibung wird eine EWMA-\brgehens- 
weise bei der PLS-Q- und der EPLS-Q-Statistik angewendet. Aus diesen Griinden wird jede \fertrauensgrenze empirisch 

40 in der vorliegenden Beschreibung bestimmt, wie dies bereits vorgeschlagen wurde [Box, 1978]. 

[0057] Die Diagnose des festgestellten, anomalen Prozessverhaltens kann durch Untersuchung der Residuen der Re- 
sponsvariablen durchgefuhrt werden. Diese Residuen konnen zu jedem Zeitpunkt in einem Balkendiagramm aufgetragen 
werden. Von einem groBen Residuum einer bestimmten Responsvariablen wird angenommen, dass es durch das anomale 
Prozessverhalten beeinflusst wird, und umgekehrt. Weiterhin fiihrte [Kourti, 1995] aus, dass ein Balkendiagramm auch 

45 durch die Residuen der Pradiktorvariablen zu jedem Zeitpunkt erzeugt werden kann. Wenn das Residuum einer bestimm- 
ten Pradiktorvariablen groBer ist, dann wird auch bei dieser Variablen angenommen, dass sie ebenfalls durch das an- 
omale Prozessverhalten beeinflusst wird. Die "GroBe" der Residuen einer bestimmten Pradiktor- oder Responsvariablen 
ist relativ in Bezug auf die Residuen anderer Pradiktor- und Responsvariablen. Ein Vergleich der momentanen Residuen 
muss ebenfalls relativ zu den Residuen der Bezugsdaten durchgefuhrt werden. Die Balkendiagramme werden auch als 

50 Fehlerbeitragsdiagramme (EC-Diagramme) bezeichnet, eines fur die Responsvariablen und eines fur die Pradiktorvaria- 
blen. Die Balkenhohen reprasentieren daher die quadrierten Residuen der Responsvariablen und die quadrierten Resi- 
duen der Pradiktorvariablen. Um diese Werte miteinander statistisch vergleichen zu konnen, muss eine Normierung 
durchgefuhrt werden. Erfolgt dies nicht, verursacht eine Responsvariable, die nicht so gut vorhergesagt werden kann wie 
andere, beispielsweise im Mittel groBere Balken in Bezug auf die anderen Responsvariablen, und umgekehrt. 

55 

4. FALLUNTERSUCHUNGEN 



4.1 Katalytische Fluidkrackeinheit 

60 [0058] Eine katalytische Fluidkrackeinheit oder FCCU ist eine okonomisch bedeutsame Einheit bei Raffiniervorgan- 
gen. Sie empfangt typischerweise mehrere verschiedene, schwere Rohstoffe von anderen Raffiniereeinheiten und krackt 
diese Strome zur Erzeugung leichterer, wertvollerer Komponenten, die schlieBlich zu Benzin und anderen Produkten ge- 
mischt werden. Die spezielle Einheit Modell IV, die von [McFarlane, 1993] beschrieben wurde, ist in Fig. 1 dargestellt. 
Der hauptsachliche Rohstoff der Einheit ist Gasol, jedoch tragen auch Strome aus schwererem Diesel und Waschol zum 

65 gesamten Rohstoffstrom bei. Frischer Rohstoff wird in einem Warmetauscher und einem Ofen vorerwarmt und dann 
dem Riser zugefuhrt, indem er mit heiBem, regeneriertem Katalysator von dem Regenerierer gemischt wird. Eine Auf- 
schlammung von dem Hauptfraktionierkolonnenboden wird ebenfalls zum Riser recycelt Der heiBe Katalysator stellt 
die Warme zur Verfugung, die fiir die endothermen Krackreaktionen benotigt wird Die gasformigen, gekrackten Pro- 
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dukte werden an die Hauptfraktionierkolonne zur TVennung geleitet Feuchtgas von der Oberseite der Hauptfraktionier- 
kolonne wird auf den Druck der Light-En ds-Anlage durch den Feuchtgaskompressor erhoht, Eine weitere Trennung 
leichter Komponenten tritt in diesem Light-Ends-Trennabschnitt auf. 

[0059] Die ausgewahlten Pradiktorvariablen fur die FCCU-Fallstudie sind in Tabelle 1 angegeben. 

TabeUe 1 



Ausgewahlte Pradiktorvariablen fur die FCCU-Fallstudie 



PRADIKTORVARIABLEN 


SIGNAL 


WaschSlf lufirate 


zeitlich konstant 


Diesel-Flussrate 


ARMA-Sequenz 


gesamter frischer Rohstoff 


ARIMA-Sequenz 


Auf schlammungs-Flussrate 


ARIMA-Sequenz 


Vorheizer-Auslass tempera tur 


ARIMA-Sequenz 


Reaktoreinstellpunkt 


konstant 


Feuchtgaskompressor- 
Saugventilposition 


abhangig vom Reaktordruck 



[0060] Alle diese Variablen gehoren zum RohstofFabschnitt der Einheit Um realistische Storbedingungen zu simulie- 
ren, wurden verschiedene Signale des typs Autoregressive Integrated Moving Average (ARIMA) diesen Variablen iiber- 
lagert, mit der Ausnahme der Dieselflussrate, die nur ein Signal des typs Autoregressive Moving Average (ARMA) 
empfing, und mit Ausnahme der Wascholflussrate und des Reaktoreinstellpunkts, die immer konstant oder gleich Null 30 
waren. 

[0061] Die Responsgruppe umfasste uberschiissigen Sauerstoff in dem Rauchgas, Konzentration von Kohlenmonoxid 
in dem Rauchgas, Riser-Temperatur, Regeneratorbetttemperatur, Regeneratorsteigrohrpegel sowie neun weitere gemes- 
sene Variablen des Systems, vergleiche [McFarlane, 1993] beziiglich einer vollstandigen Liste gemessener Variablen fur 
das FCCU-System. 35 
[0062] Um die PLS- und EPLS- Algorithmen zu testen, wurde der FCCU-Simulator vergroBert, so dass er verschiedene 
vorprogrammierte Fehler enthielt, die auf Befehl eingesetzt werden konnten. Der erste war eine stufenweise Anderung 
auf den Koksablagerungsfaktor des Rohstoffs, wodurch ein Pfropf aus einem schwerer als normalen Rohstoff simuliert 
wurde, der in die Einheit hineingelangt. Der zweite simulierte eine Unterbrechung des Flusses des regenerierten Kataly- 
sators zwischen dem Regenerierer und dem Riser, was typischerweise durch teilweise oder vollstandige Verstopfung von 40 
Dampfinjektoren hervorgerufen wird, die sich in dieser Leitung befinden. 

[0063] Beim ersten Durchlauf war kein fortgeschrittenes Regelsystem vorhanden, nur Regelsteuerungen fur den Reak- 
tor, Luftkompressorflussraten und den Reaktordruck. Bei dieser Regelkonfiguration war keine Riickkopplung zwischen 
Respons- und Pradiktorvariablen bei den Laufen 1 und 2 vorhanden, da die Feuchtgaskompressor-Saugventilposition 
weggelassen wurde. 45 
[0064] Fig, 3 zeigt die PLS-T 2 - und Q-Uberwachungsdiagramme sowie die entsprechenden EPLS-Diagramme uber ei- 
nen Zeitraum von annahernd 1500 Stunden im Normalbetrieb. In alien Figuren, in denen T 2 - und Q-Plots gezeigt sind, 
gibt die obere durchgezogene Linie die Vertrauensgrenze von 99% fur die jeweilige aufgetragene Statistik an, wahrend 
die untere gestrichelte Linie die Vertrauensgrenze von 95% angibt. Weiterhin ist die Ordinate jedes T 2 - und Q-Plots im 
Zehnerlogarithmus dargestellt Die Abtastperiode wurde auf 30 Minuten eingestellt. 50 
[0065] In Fig. 4 sind die Responses der T 2 - und Q-Statistiken von PLS und EPLS fur den ersten Fehler gezeigt, der bei 
annahernd 190,5 Stunden auftritt. Da dieser Fehler eine Anderung der Zusammensetzung des Rohstoffes simuliert, nam- 
lich einen Pfropf aus schwererem Rohstoff, wirkt er sich sofort in dem Riser aus, und danach in anderen Teilen der Ein- 
heit, die durch eine Anderung des Riser-Zustands beeinflusst werden, nach 191 Stunden. Allerdings gibt es keinen direk- 
ten, mechanischen Weg zuriick zu einem Teil des Rohstoffsystems, und daher wird keine der Pradiktorvariablen, wie sie 55 
fiir die Laufe 1 und 2 definiert wurden, direkt beeinflusst. Ebenso wenig werden sie durch Riickkopplung von den Re- 
sponsvariablen beeinflusst, da kein fortgeschrittenes Regelsystem vorhanden ist, um eine derartige Riickkopplung zur 
Verfugung zu stellen. Daher ergibt sich aus den PLS-T 2 - und Q-Diagrammen iiberhaupt keine Anzeige, dass der Fehler 
einen anomalen Zustand hervorgerufen hat. 

[0066] Im Gegensatz hierzu identifizieren die EPLS-T 2 - und Q-Statistiken, die in Fig. 4 aufgetragen wurden, deutlich 60 
den anomalen Zustand an der Vertrauensgrenze von 99%. Das EC-Diagramm, entsprechend der Zeit, zu welcher dieses 
Ereignis deutlich wird (nach 191 Stunden), ist in Fig. 5 gezeigt Die Variablen 12 und 13, namlich uberschussiger Sau- 
erstofT in dem Rauchgas bzw. Konzentration von Kohlenmonoxid in dem Rauchgas, tragen deutlich zu dem Ereignis bei. 
Dies macht physikalisch Sinn, da ein Pfropf aus einem schwereren Rohstoff zu einem rapiden Anstieg der Menge an 
Koks fuhrt, die auf dem Katalysator in dem Riser abgelagert wird, und zum Regenerierer transportiert wird, mit einer di- 65 
rekten Auswirkung auf den Sauerstoffverbrauch und die Erzeugung von Kohlenmonoxid. Das Beitragsdiagramm ergibt 
nicht direkt die potentielle Quelle des Fehlers, sondern stellt einem erfahrenen Anlagenbetreiber Information zur \ferfu- 
gung, die dazu hilft, potentielle Ursachen einzugrenzen. Im Gegensatz zu dem EC-Diagramm fur die Responsvariablen 
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zeigt das EC-Diagramm fur die Pradiktorvariablen keinen groBen Beitrag ftir iigendeine Variable. 
[0067] Bei dem zweiten Lauf wurde der Fehler des regenerierten Katalysators nach 322 Stunden eingesetzt Da die 
Pradiktorvariablen alle aus dem RohstofFabschnitt der Einheit stammen, ergibt wiederum ein Fehler oder eine Storung, 
die entweder im Reaktor oder dem Regenerierer auftreten, oder in den Verbindungskatalysatorleitungen, keinen mecha- 

5 nistischen Weg zuriick zu diesen Variablen. In diesem Fall beeinflusst der Fehler nur Responsvariablen und stellen her- 
kommliche PLS-T 2 - und Q-Diagramme das Ereignis nicht fest. Dies wird aus den oberen beiden Plots von Fig. 6 deut- 
lich. Die EFLS-T 2 - und Q-Diagramme stellen jedoch den anomalen Zustand deuthch nach 322,5 Stunden fest Das ent- 
sprechende EC-Diagramm fur die Responsvariablen, das in Fig, 7 gezeigt ist, zeigt an, dass iiberschussiger Sauerstoff in 
dem Rauchgas und der Standrohrpegel in dem Regenerierer signifikante Beitrage liefern. Dies lasst sich einfach erklaren, 

10 da jede Anderung des Russes des regenerierten Katalysators das Materialgleichgewicht in dem Standrohr beeinflusst, 
und daher dessen PegeL Eine Anderung des Russes des regenerierten Katalysators beeinflusst auch das \ferhaltnis von 
Katalysator zu Rohstoffin dem Riser, was zu einer Anderung der Menge an Koks fuhrt, die auf verbrauchtem Katalysa- 
tor abgelagert wird, und dann zu einer Anderung des Pegels des Sauerstoflverbrauchs in dem Regenerierer. Es wird dar- 
auf hingewiesen, dass das EC-Diagramm fur die Pradiktorvariablen keinen anomal groBen Beitrag fur irgendeine Varia- 

15 ble zeigt, da keine Regelungsriickkopplung vorhanden ist. 

[0068] Fur die Laufe 3 und 4 wurde die Position des Feuchtgaskompressorsaugventils als Pradiktorvariable mit be- 
riicksichtigt. Die Auswirkung irgendeiner Storung oder irgendeines Fehlers, welche den Reaktordruck beeinflussen, wird 
daher an die Pradiktorvariable ubertragen, die durch die Ruckkopplungswirkung der Reaktordruckregelung festgelegt 
wird. In diesem Fall wiirde man sowohl bei PLS als auch bei EPLS erwarten, dass sie einen anomalen Zustand feststel- 

20 len, der sich bei dieser Art von Fehler ergibt, und dies zeigt sich deuthch in Fig. 8 fur den ersten Fehler, und in Fig. 9 fur 
den zweiten Fehler. 

[0069] Der Einsatz der Vorgehensweise I bei der FCCU-Falluntersuchung hat gezeigt, dass die PLS-Uberwachungs- 
diagramme nur ein anomales Prozessverhalten feststellen konnen, wenn eine Steuerungsriickkopplung bei den Pradik- 
torvariablen vorhanden ist. Dies fuhrt weiter dazu, dass die EC-Diagramme fur die Pradiktorvariablen keinen anomal 
25 groBen Beitrag zu irgendeiner Variablen zeigen. Im Gegensatz hierzu stellten die EPLS-Diagramme beide Fehler fest. 
Beim Vorhandensein einer Steuerungsriickkopplung sind jedoch auch die PLS-Diagramme empfindlich. 

4.2 FlieBbettreaktor 

30 [0070] Dieser industrielle Prozess erzeugt zwei Losungsmittelchemikalien, die mit F und G bezeichnet sind, und be- 
steht aus mehreren Betriebseinheiten. Die Hauptelemente dieser Anlage sind kontinuierlich arbeitende Einheiten, in de- 
nen die chemische Reaktion durchgefuhrt wird. Diese Einheiten sind funf parallel arbeitende FlieBbettreaktoren, in de- 
nen jeweils F und G durch komplizierte, exothermische, chemische Reaktionen erzeugt werden. Diese Reaktoren werden 
mit zwei unterschiedlichen Stromen aus fUnf verschiedenen Reaktanten versorgt. Fig. 2 zeigt schematisch einen Reaktor 

35 und seine benachbarten Einheiten. 

[0071] Der erste Strom besteht aus dem Reaktanden A, und der zweite Strom aus den Reaktanten B, C, D und E. A und 
B sind Molekiile des typs X 2 , C ist eine Saure, D sind Molekule, die in stromaufwartigen Einheiten erzeugt werden, und 
E sind recycelte Erzeugnisse der Anlage. D und E werden durch einen stromaufwartigen Verdampfer verdampft, bevor 
sie in den Reaktor als Teil des zweiten Stroms hineingelangen. 

40 [0072] SchheBUch wird nach Verlassen des Reaktors die Trennung von F und G durch stromabwartige Destillierein- 
heiten erreicht. 

[0073] Der Reaktor besteht aus einem groBen Mantel und einer Anzahl vertikal angeordneter Rohre, in denen die che- 
mische Reaktion durchgefuhrt wird, unterstiitzt durch einen FlieBbettkatalysator Ein Thermoelement ist am Boden jedes 
Rohrs angebracht, um die Temperatur des FlieBbettes zu messen. Zum Abfuhren der Warme der exothermen Reaktion 
45 wird Ol um die Rohre herum umgewalzt. Das Verhaltnis von F : G wird regelmaBig in einem Labor analysiert Auf der 
Grundlage dieser Analyse wird das Verhaltnis von F : G durch die ReaktorrohstofFzufuhrraten eingestellt. Damit der Ka- 
talysator standig im FlieBbett vorhanden ist, wird die FlieBbettgeschwindigkeit dadurch konstant gehalten, dass der Re- 
aktordruck in Bezug auf die Gesamtflussrate eingestellt wird. 

[0074] Die chemische Reaktion wird durch nicht gemessene Storungen und Anderungen des FlieBbettes des Katalysa- 
50 tors beeinflusst. Die am haufigsten beobachtete, nicht gemessene Storung wird durch Umschlagen des Drucks des 
Dampfflusses hervorgerufen, der vom Verdampfer benotigt wird. Ungemessene Storungen konnen auch durch das Kiihl- 
mittel (Ol) hervorgerufen werden, das von einer getrennten Einheit zur Verfugung gestellt wird. Infolge des Regelsche- 
mas der Verdampfereinheit fuhren die Druckumschlage des Dampfflusses zu einer groBeren oder kleineren Russrate des 
zweiten Stroms, der in den Reaktor eintritt FlieBbettprobleme treten auf, wenn die Katalysatorverteilung uber ein Rohr 
55 am Boden des Rohrs erheblich groBer ist. Dies bringt es mit sich, dass die chemische Reaktion am Boden des Rohrs ver- 
starkt wird, was zu einer signifikanten Erhohung der Rohrtemperatur fuhrt. 

[0075] Wahrend eines Zeitraums von mehreren Wochen wurden normale Betriebsdaten sowie Daten, welche Prozes- 
sanomalitaten enthielten, fur einen bestimmten Reaktor erhalten. Die Datengruppe zum Aufhehmen des normalen Pro- 
zessbetriebs (Bezugsdaten) musste sorgfaltig ausgewahlt werden. Es musste sichergestellt werden, dass die Bezugsdaten 
60 keine Storungen einfangen, wie sie voranstehend beschrieben wurden, oder FlieBbettprobleme eines oder mehrerer 
Rohre. Weiterhin konnte, wenn die Menge an Bezugsdaten zu gering ware, eine normale Variation, die wahrend der che- 
mischen Reaktion auftritt, nicht in der Beschreibung enthalten sein. Jede Datengruppe beschreibt den Prozess im Dauer- 
zustandsbetrieb. 

[0076] Zum Identifizieren eines Dauerzustands-FLS-Modells mussten Pradiktor- und Responsvariablen ausgewahlt 
65 werden. Die Pradiktorvariablen sind die Flussrate der Reaktanten A, B, D und E, der Dampffluss zum Verdampfer, und 
ein zusatzlicher Strom, der zur Verringerung des Drucks in dem Verdampfer erforderlich ist. Als Responsvariable wurde 
die Temperatur jedes Rohrs ausgewahlt 

[0077] Eine Voruntersuchung der Daten ergab, dass die Rohrtemperaturen stark korreliert sind. Weiterhin ist eine Kor- 
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relation auch zwischen den Pradiktorvariablen vorhanden. Die Bestimmung der Anzahl an LVs, die beibehalten werden 
sollen, ergab jedoch, dass samtliche sechs LVs significant zur Vorhersage der Responsmatrix beitragen. Die Auswahl, 
wieviele LVs beibehalten werden sollen, wurde dadurch durchgefuhrt, dass eine Kreuzuberprufung so durchgefuhrt 
wurde, dass jeweils eine GroBe weggelassen wurde. Dieser Fall stellt daher ein Beispiel fur den zweiten Nachteil dar, der 
auftritt, wenn die Prozesszustandsuberwachung mit der \brgehensweise I durchgefuhrt wird. 5 
[0078] Obwohl jeder Ergebnisvektor der Pradiktormatrix beibehalten werden muss, verringert PLS die Anzahl an Pro- 
zessvariablen betrachtlich. Nachdem die Idenufizierungsprozedur beendet wurde, wurden die sechs verallgemeinerten t- 
Ergebnisse und die verallgemeinerten Residuenergebnisse fur die Bezugsdaten gemaB Gleichung (9) berechnet. Hieran 
schloss sich die Berechnung der entsprechenden Werte der PLS-T 2 -, EPLS-T 2 - und Q-Statistik und der Schwellenwerte 
fiir die zugehorigen tJberwachungsdiagramme an. Es wird darauf hingewiesen, dass die PLS-Q-Statistik nicht bestimmt 10 
werden kann. 

[0079] Fig. 1 2 zeigt die Uberwachungsdiagramme fur die PLS- und die EPLS-\brgehensweise fur die Bezugsdaten. Es 
wird darauf hingewiesen, dass die Werte jeder Statistik im logarithmischen MaBstab dargestellt sind. Die Diagramme der 
PLS-T 2 - und der EPLS-T 2 -Statistik zeigen eine natiirliche Variation des Prozesses, beispielsweise infolge von Anderun- 
gen der Zufuhr. Weiterhin zeigt das Diagramm der EPLS-T^-Statistik den Einfluss gemeinsamer Variationen (nicht ge- 15 
messener Storungen), die das Modell nicht beschreiben kann. 

[0080] Das erste beobachtete, anomale Prozessverhalten reprasentiert eine groBe, nicht gemessene Stoning infolge des 
Ab falls des Dampfdrucks. Die sich ergebenden Oberwachungsdiagramme sind in Fig, 13 dargestellt. Obwohl die 
Dampfrate konstant bleibt, andert sich das Enthalpiegleichgewicht innerhalb des Verdampfers und beeinflusst die Zu- 
sammensetzung der Reaktanten D und E in dem zweiten Strom zum Reaktor. Die nicht gemessene Storung trat nach etwa 20 
1300 Minuten in der aufgezeichneten Datengruppe auf. Die EPLS-Q-Statistik stellt diese nicht gemessene Storung un- 
mittelbar danach fest, da die Zusammensetzung des zweiten Stroms deutlich die Reaktionsbedingungen beeinflusst, wel- 
che das PLS-Modell nicht beschreiben kann. Die nicht gemessene Storung wird nicht von der PLS-T*- und der EPLS -T 2 - 
Statistik aufgenommen, da die Prozessvariation dieses Ereignisses nicht die Variationen in den Bezugsdaten uberschrei- 
tet. Die entsprechenden EC-Diagramme fur 1500, 1501 und 1502 Minuten, vergleiche Fig. 4, zeigen, dass bei etwa der 25 
Halfte der Rohre die Temperatur nicht exakt in Bezug auf die Bezugsdaten vorhergesagt werden kann. Die nicht gemes- 
sene Storung beeinflusst deutlich die Reaktionsbedingungen in samtlichen Rohren, was dadurch bestatigt werden konnte, 
dass hintereinander die drei EC-Diagramme untersucht wurden. Die Diagnose dieses anomalen Prozessverhaltens ergibt 
sich fur einen erfahrenen Anlagenbetreiber, der aus der bereitgestellten Information zuriick auf den Abfall des Dampf- 
drucks schlieBen kann. 30 
[0081] Das zweite anomale Prozessverhalten beschreibt ein FlieBbettproblem in einem der Rohre. Es gibt einige Ma- 
nipulationen, die ein Anlagenbetreiber durchfuhren kann, um das FlieBbett zu verbessem und daher die Rohrtemperatur 
zuriick auf ihren normalen Betriebswert zu bringen. Die erste Temperaturerhohung wurde jedoch von dem Anlagenbe- 
treiber nicht bemerkt. Als die zweite Temperaturerhohung festgestellt wurde, wurde der Versuch untemommen, die Tem- 
peratur zuriick zum Normalbetrieb zu bringen. Fig. 15 zeigt die entsprechenden PLS-T 2 - und EPLS-Uberwachungsdia- 35 
gramme. Beide EPLS-Statistiken stellen in beiden Fallen fest, dass die Rohrtemperatur anomal hoch ist. Im Gegensatz 
hierzu lost die PLS-T 2 -Statistik nur einen Alarm an der Vertrauensgrenze von 99% nach dem zweiten Temperaturanstieg 
aus. Auch die PLS-T 2 -Statistik uberschreitet zumindest die Vertrauensgrenze von 95%, und zeigt daher ein anomales 
Prozessverhalten an. Die Empfindlichkeit des PLS-T^Diagramms ist deswegen vorhanden, da die Ruckkopplung der 
Regelschleife fur die FlieBbettgeschwindigkeit infolge des anomalen Rohrverhaltens vorhanden ist. Die Manipulation 40 
des Prozessbetreibers lasst sich aus dem deutlichen Kurvenknick in jedem Uberwachungsdiagramm etwa bei 900 Minu- 
ten in Fig. 15 erkennen. Nachdem dieser Versuch fehlschlug, schaltete der Anlagenbetreiber schlieBlich das Rohr ab. Die 
Oberwachungsdiagramme der EPLS-T 2 - und der Q-Statistik entsprechen dem Abschalten durch Stillegen. Fig. 16 zeigt 
die EC-Diagramme nach 436 Minuten, 880 Minuten, 905 Minuten und 910 Minuten. Entsprechend den EC-Diagrammen 
konnte der Anlagenbetreiber mehr als 450 Minuten fruher MaBnahmen eingeleitet haben, um den Betrieb dieses Rohrs 45 
aufrechtzuerhalten. 

[0082] Die Anwendung der Vorgehensweise I bei diesem industriellen Beispiel zeigte, dass PLS-Uberwachungsdia- 
gramme unempfindlich (zumindest beim ersten Beispiel) in der Hinsicht sind, ein anomales Verhalten des Prozesses fest- 
zustellen. Bei dem zweiten Beispiel konnte jedoch PLS das anomale Rohrverhalten feststellen (zumindest bei der Ver- 
trauensgrenze von 95%). Im Gegensatz hierzu konnten die EPLS -Uberwachungsdiagramme deutlich jede Prozessano- 50 
malitat feststellen (an der Vertrauensgrenze von 99%). Weiterhin musste jeder Wert fur LV beibehalten werden, infolge 
seiner Signifikanz zum Vorhersagen der Responsvariablen. Daher konnte nur ein Uberwachungsdiagramm bei der \for- 
gehensweise I erhalten werden. 

[0083] Die zweite Vorgehensweise (Vorgehensweise II) wurde zu 16 Streuplots fuhren und sechs Oberwachungsdia- 
grammen, von denen keines eine Variation der Responsvariablen beschreibt. Im Gegensatz hierzu benotigt EPLS nur 55 
zwei Uberwachungsdiagramme und beschreibt die EPLS-T 2 -Statistik eine Variation der Responsvariablen, vergleiche 
Gleichung 9. 

5. SCHLUSSFOLGERUNGEN 

60 

[0084] Bei der vorliegenden Beschreibung werden die herkornmlichen PLS-Vorgehensweisen fur die Zustandsuberwa- 
chung kontinuierlicher, industrieller Prozesse, wie sie in [MacGregor, 1991, 1995], [Kresta, 1991] und [Wise, 1996] be- 
schrieben werden, erneut uberpnift, und werden Problembereiche hervorgehoben. Aus dieser Untersuchung ergibt sich, 
dass die herkornmlichen PLS-Oberwachungsdiagramme entweder unempfindlich oder schwierig zu analysieren sein 
konnen, falls sich der Prozess anomal verhalt. Diese Beschreibung stellt eine Ausweitung des Standard-FLS-Algorith- 65 
mus zur Verfugung, bezeichnet als EPLS, die zu der Definition neuer PLS-Ergebnisse fuhrt, die als verallgemeinerte Er- 
gebnisse bezeichnet werden. Ahnlich wie bei herkonunlichen PLS-Vorgehensweisen, konnen Statistiken auf der Grund- 
lage der verallgemeinerten Ergebnisse von EPLS definiert werden, die in Abhangigkeit von der Zeit auf Uberwachungs- 
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diagrammen aufgetragen werden konnen. Diese Oberwachungsdiagramme beschreiben die Gesamtvariation der Pradik- 
tor- und Responsvariablen (EPLS-T^Diagramm) und ihrer Residuen (EPLS-Q-Diagramm). 

[0085] Eine theoretische Analyse der Oberwachungsdiagrarnrne, die aus den verallgemeinerten Ergebnissen von EPLS 
abgeleitet werden, und herkommlicher PLS-Vorgehensweisen ergibt Folgendes: 

5 

1. Bei der Vorgehensweise I wind, wenn ein anomales Verhalten Responsvariablen beeinflusst, die nicht Teil einer 
geschlossenen Regelschleife sind, diese Situation moglicherweise nicht festgestellt Bei EPLS wird die anomale Va- 
riation der Responsvariablen deutlich. 

2. Falls die Responsvariablen stark korreliert sind, jedoch die Pradiktorvariablen nicht, erzeugt die \forgehensweise 
10 I nur ein Diagramm, das T 2 -Diagramm. Bei EPLS bleiben beide Diagramme relevant, unabhangig von der Anzahl 

an latenten Variablen, die beibehalten werden. 

3. Die zweite Vorgehensweise (Vorgehensweise H) definiert eine Anzahl an Diagrammen in Abhangigkeit von der 
Anzahl an latenten Variablen, die beibehalten werden. Im Gegensatz hierzu benotigt EPLS nur zwei Diagramme, 
unabhangig von der Auswahl ausgewahlter, latenter Variablen. 

15 4. Der Einsatz des VorfiersagefeMerquadratdiagramms zusammen mit den beiden ttberwachungsdiagrammen der 

Vbrgehensweise I kann die voranstehend geschilderten Schwierigkeiten uberwinden. Allerdings fiihrt dies zu zu- 
mindest drei Uberwachungsdiagrammen, und ist die Variation der Responsvariablen in keinem dieser Diagramme 
vorhanden. Im Gegensatz hierzu benotigt EPLS nur zwei Uberwachungsdiagramme, und wind die Variation der Re- 
sponsvariablen nicht in dem EPLS-I^-Diagramm akkumuliert 

20 

[0086] Die vorliegende Beschreibung stellt auch zwei Einsatzuntersuchungen zur Verfugung, um die voranstehenden, 
theoretisch abgeleiteten Ergebnisse zu bestatigen. Diese Anwendungen betreffen die Simulation einer katalytischen 
Ruidkrackeinheit (FCCU) und einen realen industriellen Prozess. Zwei anomale Situationen sind in beiden Falluntersu- 
chungen vorhanden, welche die Auswirkung einer nicht gemessenen Stoning und eine "interne" Anderung des Prozess- 
25 verhaltens beschreiben. 

[0087] Die Ergebnisse der FCCU-Falluntersuchung zeigen deutlich, dass eine Regelungsriickkopplung fur die \fcrge- 
hensweise I unbedingt erforderlich ist, um ein robustes und empfindliches Zustandsuberwachungswerkzeug zur Verfu- 
gung zu stellen. Falls dies nicht garantiert ist, kann es bei dieser \forgehensweise auftreten, dass anomale Situationen 
nicht festgestellt werden. Im Gegensatz hierzu stellt die EPLS-\brgehensweise robuste und empfindliche Uberwa- 

30 chungsdiagramme zur Verfugung, unabhangig von dem Vorhandensein einer Regelungsriickkopplung. 

[0088] Dies kann auch durch die Untersuchung des industriellen Falls bestatigt werden. Die erste anomale Situation 
wird nicht von der PLS-T 2 -Statistik festgestellt, da der Umschlag des Dampfdrucks nicht die Pradiktorvariablen stark ge- 
nug beeinflusst, um festgestellt zu werden. Es wird allerdings der Zustand der chemischen Reaktion innerhalb der Rohre 
beeinflusst und daher das Verhalten der Responsvariablen. Die zweite Situation beschreibt ein anomales Verhalten eines 

35 der Rohre. In dieser Situation beeinflusst die Regelungsriickkopplung die Pradiktorvariablen infolge des anomalen \fer- 
haltens. 

[0089] Die Untersuchung des industriellen Falls zeigt weiterhin, dass samtliche latenten \fcriablen signiflkant zur Vor- 
hersage der Responsvariablen beitragen. Daher ist nur ein Uberwachungsdiagramm (das T 2 -Diagramm) fur die Vorge- 
hensweise I verfugbar. Weiterhin wurde die Vorgehensweise II zu insgesamt 16 Streuplots und 6 Uberwachungsdiagram- 
40 men fuhren. Bei dieser Anzahl an relevanten Diagrammen ist es entsprechend muhsam, ein anomales Verhalten festzu- 
stellen, im Unterschied zur Beobachtung der beiden EPI^-Uberwachungsdiagrarnme. 

[0090] Die weitere Forschung in Bezug auf die verallgemeinerten Ergebnisse von EPLS konzentriert sich auf Anwen- 
dungen, welche dynamische Prozessmodelle verwenden, wie sie beispielsweise fur eine Modellvorhersageregelung er- 
forderlich sind. Es soil untersucht werden, ob dieses dynamische Prozessmodell auch die Grundlage fur die Prozesszu- 
45 standsiiberwachung bereitstellt, im Unterschied zu der Dauerzustandsuntersuchung, die bislang vorgenommen wurde. 
Dariiber hinaus wird der Einsatz der verallgemeinerten EPLS-Ergebnisse fur die Uberwachung von postenweise ab- 
laufenden Prozessen der Gegenstand zuktinftiger Uberlegungen sein. Die PLS-\forgehensweise zur Uberwachung po- 
stenweise ablaufender Prozesse, die in [Nomikos, 1994 und 1995] diskutiert wurde, wird daher die Grundlage fur die Un- 
terscheidung eines "guten" Postens von einem "schlechten" Posten sein. 

50 

ANHANGE 



A.l Theoretische Analyse des PLS-Algorithmus 

55 [0091] Der PLS-Identifizierungsalgorithmus beruht auf der Bestimmung jedes Paars von Komponentenmatrizen 
und Yk (siehe Gleichung 1), durch eine iterative Prozedur. Nachdem der k-te Iterationsschritt durchgefuhrt wurde, wer- 
den die berechneten Komponentenmatrizen von der Pradiktormatrix bzw. der Responsmatrix subtrahiert, vor derBerech- 
nung des (k+l)-ten Iterationsschrittes. Die Subtraktion der Komponentenmatrizen wird auch als Druckablassprozedur 
bezeichnet und verlauft folgendermaBen: 

60 _ ^ 

Y, -Y w -fi,qr * (A1) 

65 [0092] Die Ergebnisvektoren t^, Cik und Belastungsvektoren p^ und werden dazu bestimmt, um den Beitrag jedes 
Paars von Komponentenmatrizen zu den Pradiktor- und Responsmatrizen zu maximieren. Dies wird dadurch erreicht, 
dass folgende Kriterien erfullt werden: 
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[0093] Losungen fur die drei Kostenfiinktioneri miissen daraufhin berechnet werden. Beginnend mit Gleichung (A2) 
werden Wk und Vk als Gewichtsvektor der Pradiktormatrix bzw. der Responsmatrix bezeichnet, und reprasentiert den 
Wert der entsprechenden Kostenfunktion. Nach [Hoskuldsson, 1988] kann einschlieBlich der Randbedingungen fiir die 
Lange der Gewichtsvektoren Gleichung (A2) folgendermaBen umgeschrieben werden: 25 

X|_|Y A _,Y J £,X > _ l w 4 = 

Y r X X r Y v - 1 v " (A5) 

[0094] Aus Gleichung (A5) wird deutlich, dass der Eigenvektor ist, der dem groBten Eigenwert der Kreuzcovari- 
anzmatrix 
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X Y<-iY Xk-i 

k-l k-l 35 

zugeordnet ist und Vk der Eigenvektor ist, der dem groBten Eigenwert der Kreuzcovarianzmatrix 

T 

Y X:<- 

k-l 40 
T 

iX f Yk-i 
it — 1 

45 

entspricht. In Gleichung (A3) ist bk der Regressionskoeffizient zwischen dem k-ten Paar von Ergebnisvektoren, namlich 
tk und Uk, und ist J q der Wert der Kostenfunktion. Die Losung von Gleichung (A3) ist die ubliche Losung der kleinsten 
Fehlerquadrate fur b k . Fiir Gleichung (A4) ist J p die zugehorige Kostenfunktion zur Bestimmung von pk, und von J q zur 
Berechnung von q k . Die Kostenfunktionen J p und J q werden durch die folgende Losung minimiert [Geladi, 1986] 

Y ; ; • (A6) 

[0095] SchlieBlich kann unter Einfluss von n latenten Variablen die Matrix der Regressionskoeffizienten 

B M_ 60 

PLS 



zwischen der Pradiktormatrix und der Responsmatrix folgendermaBen (Lindgren, 1993] berechnet werden: 

Y D = X 0 .B{& + E„ 65 

Bja-w.fcrw.j-^b.wr (A7) 
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wobei W n , P n und Qn Matrizen sind, in denen die n Gewichtsvektoren w k und Belastungsvektoren p* und als Spalten 
gespeichert sind. Gemafi Gleichung (A2) wind der Gewichtsvektor flir die Pradiktormatrix, Wk, mit der Druckablass-Pra- 
diktormatrix X k _j multipliziert, um den Ergebnisvektor t k zu bestimmen. Allerdings gab [Lindgren, 1993] an, dass der Er- 
gebnisvektor ^ auch direkt aus der urspriinglichen Pradiktormatrix Xo folgendermafien berechnet werden kann: 

.o T. = X„[r, ... r.]=X„R.;R. = W.[p;\V.I' 

A.2 Fiir Gleichung (7) erforderlicher Beweis 

15 [0096] Die Vorhersage der Responsmatrix auf der Grundlage von n latenten Variablen, die beibehalten werden, ist fol- 
gendermaBen: 

Y B =T 0 diag{b o )Q:. (A9) 

20 [0097] Die t-Ergebnisse konnen direkt aus der urspriinglichen Pradiktormatrix berechnet werden, vergleiche Glei- 
chung (A8): 

T n = X a R n . (A10) 

25 [0098] Die PLS-Regressionsmatrix kann daher folgendermafien bestimmt werden: 

*&=Kdiag{bM. (All) 

30 [0099] SchlieBlich stellt die vorherige Multiplication mit PJ die erfoiderliche Gleichheit zur Verfiigung (siehe Anhang 
3): 

p:B^=^{b„}Qi. (Al2) 
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A. 3 Fiir Gleichung (9) erforderlicher Beweis 
[0100] Die Definition der Matrix R„, vergleiche Gleichung (A8), kann dazu verwendet werden, zu beweisen, dass gilt: 



wobei 3nxo eine Einheitsmatrix n bezeichnet. Die Elemente des Matrixprodukts sind folgendermafien definiert: 
45 Prr, = p. r ^-w 1 p^«tl-w i _ l p^ I ]wr / ;lSi,jjSn (AI3) 

[0101] Wenn i grofier ist als j, ist Gleichung (A 13) gleich Null, da Skalarprodukte zwischen pj T und w m auftreten, mit 
1 <m<j, was jeweils gleich Null ist [Hoskuldsson, 1988]. Gemafi Gleichung (A 13) konnen, falls i kleiner als j ist, die 
Faktoren des Matrixproduktes bis zum i-ten Faktor reduziert werden, was fiihrt zu: 

P>. = frf - pfwj»r>--[3 - w,.^;^, = fiidi „ , 

[0102] Gleichung (A 14) zeigt, dass der transponierte Vektor in dem Term auf der rechten Seite gleich Null ist, und da- 
her das gesamte Produkt gleich Null ist. Weiterhin fuhrt die Tatsache, dass pi T Wi = 1 ist, auch dazu, dass pi T ri = 1 ist Zu- 
sammenfassend wird aus den Gleichungen (A13) und (A14) deutlich, dass das Matrixprodukt Ro T Pn gleich einer Ein- 
heitsmatrix n • n ist. 

A.4 Orthogonalitat beider verallgemeinerter Ergebnisse 



50 



55 



60 



[0103] Um zu beweisen, dass beide verallgemeinerten Ergebnisse nicht orthogonal sind, werden die Kovarianzmatri- 
zen fur beide Ergebnisse, einschliefilich von n LVs, untersucht. Die Kovarianzmatrix der verallgemeinerten t-Ergebnisse 
65 ist in Gleichung (A15) angegeben: 
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R: 



[0104] In Gleichung (A 15) reprasentieren 

s (») ,W s .CO , s<">. sounds W 
JT XX FE FF 

B ft * a * ■ 



(A16) 
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Kovarianzmatrix der verallgemeinerten t-Ergebnisse, die Standard- t-Ergebnisse, die verallgemeinerten restlichen Ergeb- 
nisse, die Pradiktorvariablen, den Vorhersagefehler der Responsvariablen, und die Residuen der Pradiktorvariablen. Wei- 
terhin bezeichnet ste* die Kreuzkovarianzmatrix der Standard-t-Ergebnisse und der verallgemeinerten restlichen Ergeb- 
nisse, und ist se£eh die Kreuzkovarianzmatrix des Vorhersagefehlers der Responsmatrix und der Residuen der Pradiktor- 20 
matrix. Es wind darauf hingewiesen, dass nur die Kovarianzmatrix s$ diagonal ist, da die Standard-t-Ergebnisse gegen- 
seitig orthogonal sind [Hoskuldsson, 1988]. Aus Gleichung (A15) kann man schlieBen, dass die Matrizen se*e* und ste* 
im AUgemeinen nicht diagonal sind, selbst unter der Annahme, dass die Spalten der Residuenmatrizen En und F n und 
ebenso die Pradiktormatrix X aus Signalen in Form weiBen Rauschens bestehen, Im theoretischen Fall, in welchem die 
Anzahl an Responsvariablen gleich der Anzahl manipulierter Variablen ist, und der Prozess entkoppelt ist, sind jedoch 25 
die Kovarianzmatrix se ( *e* und die Kreuzkovarianzmatrix STE*vom Diagonaltyp. Daher kann die Kovarianzmatrix st*t* 
nicht vom Diagonaltyp sein, wenn eine Anzahl n an LVs zuriickbehalten wird. Wenn alle M (siehe Gleichung 1) LVs bei- 
behalten werden, dann reduziert sich Gleichung (A 15) auf: 



30 



[0105] Auf der Grundlage der Annahme, dass En aus Signalen in Form weiBen Rauschens besteht, ist die Kovarianz- 
matrix Se*e* offensichtlich vom Diagonaltyp. Infolge der vorherigen und nachfolgenden Multiplikation mit im AUge- 35 
meinen nicht-diagonalen Matrizen ergibt sich jedoch das Ergebnis, dass se*e* und daher st*t* nur symmetrisch sind. Die 
Interpretation der Gleichung (A15) und (A16) hat gezeigt, dass weder die Kovarianzmatrix se ( *e* noch s-pt* diagonal ist. 
Daher sind beide verallgemeinerten Ergebnistypen im AUgemeinen nicht orthogonal. 
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20 Patentanspruche 

1. Verfahren zum Konstruieren/Konfigurieren eines mehrdimensionalen Monitors fur statistische Prozesse durch 
eine Vorgehensweise mit partiellen kleinsten Fehlerquadraten, mit einer Konstruktion der Prozesspradiktor- und 
Responsmatrizen aus Bezugsdaten, wobei die Pradiktormatrix aus Signalen der manipulierten und gemessenen Sto- 

25 rungs- oder Ursachenvariablen des Prozesses (Pradiktorvariablen) besteht, und die Responsmatrix aus den gesteu- 

erten oder verursachten Variablen des Prozesses (Responsvariablen) besteht, Zerlegen der Pradiktor- und Respons- 
matrizen in Komponentenmatrizen des Rangs Eins, wobei jede der Matrizen aus einem Vektorprodukt besteht, bei 
welchem ein Vektor (der Ergebnisvektor) die Variation beschreibt, und der andere (der Belastungsvektor) den Bei- 
trag des Ergebnisvektors zu der Pradiktor- oder Responsmatrix, wobei die Zerlegung durch die Erzeugung einer pa- 

30 rametrischen Regressionsmatrix auf der Grundlage von Iterationen der Zerlegung der Pradiktor- und Responsmatri- 

zen durchgefuhrt wind, gekennzeichnet durch die Erzeugung eines ersten verallgemeinerten Ergebnisvektors, der 
irgendeine signifikante Variation des Prozesses einschlieBlich Variationen der Pradiktor- und Responsvariablen be- 
schreibt, und eines zweiten verallgemeinerten Ergebnisvektors, der den Vorhersagefehler des Modells der partiellen 
kleinsten Fehlerquadrate und Residuen der Pradiktormatrix reprasentiert. 

35 2. Verfahren zum Konstruieren/Konfigurieren eines mehrdimensionalen Prozessmonitors nach Anspruch 1, bei 

welchem die verallgemeinerten Ergebnisse durch Konstruktion einer vergroBerten Matrix berechnet werden, die 
hier mit Z bezeichnet ist, und folgende Form aufweist: 



40 



45 



Z=[YiX], 

wobei X die Pradiktormatrix und Y die Responsmatrix ist, und eine Ergebnismatrix T n = T n -En* konstruiert wird, 
bei welcher T n * und En* folgende allgemeine Form aufweisen: 



t; = [Y:x)[u-;r J ] ! R. 
e; = K:fJ[b-;31x 



die Spalten der Matrix T n * die verallgemeinerten t-Ergebnisse zur Verfugung stellen, und die Spalten der Matrix 
50 En* die verallgemeinerten restlichen Ergebnisse, wobei eine Einheitsmatrix M x M bezeichnet, und 

»i-> 

PLS 

55 die PLS-Regressionsmatrix ist. 

3. Mehrdimensionaler Monitor fur statistische Prozesse, der durch das Verfahren von Anspruch 1 oder Anspruch 2 
konstruiert/konfiguriert wurde, und so ausgebildet ist, dass er ein anomales Prozessverhalten durch Analyse der Re- 
siduen der Responsvariablen identifiziert. 

4. Verfahren zur Oberwachung eines Prozesses, welches das Konfigurieren eines mehrdimensionalen Monitors fur 
60 statistische Prozesse durch das Verfahren nach Anspruch 1 oder Anspruch 2 umfasst, und Identifizieren eines an- 

omalen Prozessverhaltens, zumindest teilweise, durch Analyse der Residuen der Responsvariablen. 
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F£HL£RB£ITRAGSOIAGRAMM FUR ZEITPUNKT 11460 min. DIE 02 - UNO C0-K0N2ENTRATION 
W DEM CIHTGASFLUSS HABEN DEN GROSSTEN VORHERSAGEFEHLER. 
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